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We model the cooling of open optical and optomechanical resonators via optical feedback in the 
Linear Quadratic Gaussian setting of stochastic control theory. We show that coherent feedback 
control schemes, in which the resonator is embedded in an interferometer to achieve all-optical 
feedback, can outperform the best possible measurement-based schemes in the quantum regime 
of low steady-state excitation number. Such performance gains are attributed to the coherent 
controller's ability to process non-commuting output field quadratures simultaneously without loss 
of fidelity, and may provide important clues for the design of coherent feedback schemes for more 
general problems of nonlinear and robust control. 
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Feedback control of classical dynamical systems plays 
a central role in modern engineering [T] but its quantum 
analogue, the notion of controlling a quantum system via 
feedback with a quantum or classical controller, is much 
less developed. Recent progress in modeling [2H5] and 
realizing [5HL3] quantum feedback underscores the need 
for systematic approaches to control design, as do the 
wide range of potential applications in quantum science 
and technology. 

While some of the most interesting problems in quan- 
tum feedback control are nonlinear [TBHIS], linear open 
quantum systems provide a logical first step towards 
more general problems. Working with linear systems, 
James, Nurdin and Petersen (20] [21] have utilized in- 
terconnection models based on quantum stochastic dif- 
ferential equations (QSDEs) [2"2"H2"5] to develop gener- 
alizations of the traditional T-L°° and Linear Quadratic 
Gaussian (LQG) control paradigms that allow for the 
possibility of coherent optical feedback with linear quan- 
tum controllers. Here we work within the quantum 
LQG framework of j^T] to study steady-state cooling 
of open quantum oscillators such as optical and op- 
tomechanical resonators subject to stationary heating 
and damping, with optical probing and feedback. We 
utilize numerical optimization together with fundamen- 
tal analytic results [H [26] to describe coherent feed- 
back control schemes that outperform the best possi- 
ble measurement-based schemes. We find more system- 
atic and quantitatively significant advantages of coherent 
feedback over measurement-based feedback than in linear 
control scenarios considered previously [TJl [21] , and in- 
terpret these performance gains in terms of the way that 
non-commuting field quadratures propagate through the 
feedback loop. 

We model quantum harmonic oscillators as cascad- 
able open quantum systems using the SLH framework 
[2"7l |2"5] and associated QSDEs. In the SLH framework, 
any open quantum system may be described as a triple 
G = (S, L, H) where S is a scattering matrix, L is a cou- 
pling vector and H is the Hamiltonian operator for the 
system's internal degrees of freedom. For a linear system 
with a vector of internal state variables x, the quantum 



stochastic differential equations (analogues to the Heisen- 
berg Equations for open quantum systems [22, 29]) take 
the following form [2"tJll2"T] : 



dxlt) = [Ax{t) +a]dt + B da{t) 
dd(t) = \Cx{t) + c]dt + Dda(t) 



(1) 



Here A, B, C, D, a and c are real and related to (S, L, H). 
The processes da(t) and dd(t) are Hermitian quantum 
stochastic processes for the inputs and outputs, respec- 
tively, defined by ddi = (dAi+dA\, (dAi — dA\)/i), where 
dA{t) is the quantum Wiener process [3U1 [3T] following 
the ltd rule dAi dA'- = Sijdt for vacuum inputs. 

As a first example of a quantum control system, take 
a noisy optical cavity, Fig. [T] This has the SLH model 



S = I, L 



fcia, vfc 2 a, yk^a 



H = Aa f a, (2) 



where A is the detuning of the cavity resonance frequency 
from that of a rotating frame. In the controller's absence, 
the cavity is driven by two vacuum inputs (mirrors k\ 
and fe) and one thermal input (mirror ^3), and the cav- 
ity relaxes to a thermal state where the mean photon 
number (a) depends on the noise power fc„. The ob- 
jective in this control problem is to minimize the effect of 



Control Objective: 

J = X 2 + P 2 k;, = 0.01 

= 4a*a + 2 
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(Noise, k„) k 3 = 0.01 = 0.01 
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FIG. 1: Possible Classical and Coherent Cavity Controllers 
for an Optical Cavity Plant System 
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the noise on the cavity's internal state - in other words, 
to minimize the photon number (at a) of the cavity. We 
accomplish this by sending output 1 through a control 
circuit and feeding the result back into input 2. This is 
an LQG feedback control problem, assuming negligible 
control cost. 

Once the full system is set up, with its associated 
A, B, C and D matrices, the covariance matrix ijjj = 
| {XiXj + XjXi) can be computed with the Lyapunov 
equation, Aa + <rA T + BFB T = 0, where F is the in- 
put noise covariance; the mean photon number can then 
be computed from a. Without a controller present, we 
find: 

h 



<N)/<N)„ 
LOi 



(at a) 



fcl + k 2 f K3 

Five possible controllers are sketched in Figure [T] The 
classical controllers (Homodyne and Heterodyne) work 
by measuring a quadrature from the cavity's output, or 
splitting the beam and measuring two different quadra- 
tures, and applying a feedback signal based on this mea- 
surement and the controller's internal state. The "trivial 
controller" works by feeding the output directly back into 
mirror 2 of the plant, perhaps with a phase shift. The re- 
maining two controllers function as coherent controllers 
with memory, where the feedback signal depends not only 
on the probe field, but also on its history. The OPO with 
squeezed input/output (I/O) is the most general system 
realizable for two inputs and one internal degree of free- 
dom [32]. 

The classical LQG optimization problem is convex, and 
the optimal controller can be expressed analytically in 
terms of the plant parameters and the LQR |26J. The 
optimal homodyne and heterodyne controllers act as clas- 
sical amplifiers. In the heterodyne case, the output da is 
related to the input da and vacuum noises daki, da k2 by: 

da x = £{da x + da kl . x ) + da k2iX 

da p = £,{da p - daui, p ) + da k2 , P (4) 

In other words, the controller measures the probe, am- 
plifies it, and then adds the signal to a vacuum channel. 
With this controller, the plant mean photon number is 



- k,, 



(3) 



(a f a) c i = 



k 2 sinh r\ + k%k n 



(5) 



ki + k 2 + k 3 + 2\/kik 2 sinh r) 

which is plotted in dotted black in Figure [2] Note that, 
while the controller does a good job in the classical 
regime, when (at a) is high, it is ineffective in the quan- 
tum regime when the photon number is < 1. 

Contrast this with the "trivial controller" , which feeds 
the output from mirror 1 directly into mirror 2. Rather 
than leaking photons separately, the two outputs then in- 
terfere constructively [TB] ■ This leads to an enhancement 
in the dissipation and a mean photon number 
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FIG. 2: Bottom: Cavity photon number as a function of noise 
strength fcjv, for various control schemes. Top: Photon num- 
ber relative to the no-control case. Smaller is better. 



which beats the best classical controller by a constant 
factor. 

A more sophisticated control design, the OPO with 
squeezed I/O, can be optimized numerically. We find a 
(locally) optimal solution which is equivalent to a two- 
mode squeezer: 



ddi = cosh(i])dai + smh(rj)da 2 
da 2 = smh(ri)da 2 + cosh(rj)da 2 

This gives a mean photon number of: 



(oto), = 

\ / 2— so 



k 2 sinh rj + k 3 k n 



(7) 



(8) 



ki + k 2 + k 3 + 2^fk^k 2 



(6) 



; + fc 3 + 2\/kxk 2 cosh j] 

which is also plotted in Figure [2j 

Notice just how similar Equations ^ and ^ are. 
Both the heterodyne controller and the linear amplifier 
reduce the cavity's photon number by amplifying the 
feedback signal, but also add noise to the system. For 
equivalent levels of amplification, the classical controller 
adds extra noise into the system from the measurement 
process. When k n and T) are large, this extra noise is 
negligible, but in the quantum regime where k n and r\ 
are < 1, this noise can play a major role in making the 
linear amplifier outperform the heterodyne controller. 

Of course, this need not be the global optimum. It 
is possible that other solutions, not found in our search 
of the parameter space, lead to better LQG performance 
than the two-mode squeezer. Nevertheless, it is signifi- 
cant that we have found a quantum controller which does 
better than the optimal classical controller for all values 
of fc„, particularly in the quantum regime where the ex- 
citation number is low. 
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Another model system to consider is the optomechan- 
ical oscillator. Optomechanical oscillators - mechanical 
springs that couple to an optical field via a cavity - have 
been a topic of tremendous recent interest in the physics 
community [331 131] • A central goal has been to find ways 
to exploit optomechanical coupling to cool the mechan- 
ical oscillator from ambient temperature to its ground 
state. 

Consider a single cavity coupled to a mirror on a 
spring. If we shine a laser into the input port, and adia- 
batically eliminate the internal mode, we obtain the SLH 
model 

S=1 2X 2, L=\ i Kx m ,^/il/Qb\, H = iltfb (9) 

where if is a constant depending on the cavity and the 
drive field. It is not hard to show that this gives the 
following input-output relations: 



dx m — Qp m QjQXyn 
dpra — Q/QPn 

da x = da x + 2Kx m dt 
dd p = da p 



2Kda r: 



(10) 



The state variable x m is imprinted on the output dd x , 
permitting the controller to be used as a "measure- 
ment" device via the x quadrature. Conversely, the p- 
quadrature input da p alters the state of the mirror, al- 
lowing the system to function as a "feedback" device. 

Coupling two such cavities together as in Fig. [3j we 
obtain our model plant system. The cooling objective 



Xm. Pm. or b, b* 




is to minimize the phonon number (b^b), where b is the 
spring's annihilation operator. As in the optical cavity 
case, this is an LQG control problem. 

Consider first the classical controller, black in Fig. [3] 
This controller only measures the da x quadrature, since 
no information is contained in da pi and applies a feed- 
back based on this input (such strategies outperform 
heterodyne-based classical control). When optimizing, 
we vary both the controller parameters and the couplings 
Ki , K 2 to the cavities in ^ , since these depend on input 
laser powers rather than on properties of the plant. Here 
we will assume K\ — —K-x = K, however, as this type of 
design could be implemented practically using two modes 
of a single optical cavity with comparable power circulat- 
ing in each. This symmetry in the control model also fa- 
cilitates straightforward interpretation of the advantage 
of coherent feedback. 

As before, the classical controller does a good job cool- 
ing the oscillator in the classical limit (N) 3> 1, but does 
a very poor job in the quantum limit; below a threshold 
value, it has no effect at all. The classical controller adds 
two sources of noise to the plant. First, the probe field 
adds measurement noise to the spring through the first 
cavity. Second, the feedback field adds a feedback noise 
of equal magnitude via the second cavity. These noises 
add up incoherently. In the low-/c n limit, they dominate 
the dynamics of the control system and the classical con- 
troller becomes ineffective. 

We also optimized two types of coherent controller - 
an optical cavity and an OPO. They both behave quali- 
tatively in a similar way, but the optical cavity is easier 
to understand. With the optical cavity controller, the 
feedback field is related to both the probe field and the 
cavity state, which changes with time: 



da = (— iA — K/2)adt + ^/KdA\ 
dAi — dA\ + \fna dt 



(11) 



Optimization shows that A = f2, so both the cavity 
and mirror resonate at the same frequency. Although 
not interpreted in terms of feedback control theory, it 
is straightforward to show that cavity-optomechanical 
cooling experiments such as [9TU1| have implemented an 
equivalent strategy in which the three optical cavities of 
Fig. [3] are collapsed into one. 



dA dA 




Coherent 
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FIG. 3: Control-system setup for the mechanical oscillator 
cooling problem with Q = 100, Q = 10000. Three potential 
controller designs and their relative performance. 



FIG. 4: Flow of x- and p-quadrature signals (blue and red, 
respectively) in the classical and coherent controllers. 
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The LQG-optimal classical controller relies on a 
Kalman filter that recursively estimates the plant state 
from the measured signal pQ. The coherent cavity con- 
troller can also be thought of as a Kalman filter, but one 
that preserves the coherence of the input signal dA\. In 
the classical controller the p-quadrature da p is essentially 
discarded after the measurement, while in the cavity con- 
troller the field retains its coherent properties and da p 
coming out is the same as da p going in. This correlates 
the noises in the measurement and feedback cavities such 
that the associated forces on the mechanical oscillator 
exactly cancel, as shown in Fig. [4] This cancellation of 
the measurement noise is what gives the coherent cavity 
controller its superior performance in the low phonon- 
numbcr regime. 

LQG control is not the first problem to benefit from 
this noise cancelation. Previous studies had shown that 
precision measurement schemes [35j . particularly in the 
context of LIGO [36], show similar improvements. 

The OPO controller is more general than the cavity 
controller and its performance is slightly better over the 
whole noise range. As with the optical cavity case, it is 
significant that we have found realistic coherent control 
schemes that perform significantly better than provably 
optimal measurement-based schemes. 

In this paper, we have studied the coherent control of 
linear quantum systems from an LQG perspective. In 
the systems studied, we have shown that there is always 
a quantum controller that does at least as well as the op- 
timal classical controller. In the quantum regime, when 
excitation number in the plant is of order unity, we have 
shown that the best quantum controller can do better 
- in some cases, significantly so. One could straightfor- 
wardly extend these results to non-quadratic cost func- 
tions in linear control systems. Indeed, some work has 
already been done on this matter, focusing on using co- 
herent feedback to maximize the squeezing in a cavity 
mode [15] . Viewing an optomechanical device as a con- 
trol system may also provide insight into minimizing the 
noise in optomechanical sensors. Finally, the fact that 
coherent controllers perform better because they utilize 
both quadratures of the input field may help guide the 
design of quantum controllers for nonlinear systems such 
as optical switches or error correcting codes. 

This work has been supported by the NSF (PHY- 
1005386), AFOSR (FA9550-11-1-0238) and DARPA- 
MTO (N66001-11-1-4106). Ryan Hamerly is supported 
by the NSF GRFP and a Stanford Graduate Fellowship. 
We thank Nikolas Tezak, Gopal Sarma, Dmitri Pavlichin 
and Orion Crisafulli for useful discussions. 



t I T-'. r- 1 

T Electronic address: hmabuchi@stanford.edu 
[1] K. Astrom and R. M. Murray, Feedback Systems: An In- 
troduction for Scientists and Engineers (Princeton Uni- 
versity Press, 2008). 
[2] V. P. Belavkin, Autom. Rem. Control 44, 178 (1983). 
[3] H. M. Wiseman and G. J. Milburn, Phys. Rev. Lett. 70, 
548 (1993). 

[4] A. C. Doherty et al, Phys. Rev. A 62, 012105 (2000). 
[5] H. M. Wiseman and G. J. Milburn, Quantum Measure- 
ment and Control (Springer, 2009). 
[6] W. P. Smith et al, Phys. Rev. Lett. 89, 133601 (2002). 
[7] M. A. Armen et al, Phys. Rev. Lett. 89, 133602 (2002). 
[8] P. Bushev et al, Phys. Rev. Lett. 96, 043003 (2006). 
[9] S. Gigan et al., Nature 444, 67 (2006). 

[10] O. Arcizet et al, Nature 444, 71 (2006). 

[11] D. Kleckner and D. Bouwmeester, Nature 444, 75 (2006). 

[12] H. Mabuchi, Phys. Rev. A 78, 032323 (2008). 

[13] G. G. Gillett et al, Phys. Rev. Lett. 104, 080503 (2010). 

[14] C. Sayrin et al, Nature 477, 73 (2011). 

[15] S. Iida et al, arXiv: 1103. 1324 (2011). 

[16] J. Kerckhoff, H. I. Nurdin, D. S. Pavlichin and 
H. Mabuchi, Phys. Rev. Lett. 105, 040502 (2010). 

[17] J. Kerckhoff, D. S. Pavlichin, H. Chalabi and H. Mabuchi, 
New J. Phys. 13, 055022 (2011). 

[18] H. Mabuchi, Appl. Phys. Lett. 98, 193109 (2011). 

[19] J. E. Gough, M. R. James, H. I. Nurdin and J. Combes, 
arXiv:1107.2976v2 (2011). 

[20] M. R. James, H. I. Nurdin and I. R. Petersen, IEEE 
Trans. Auto. Control 53, 1787 (2008). 

[21] H. I. Nurdin, M. R. James and I. R. Petersen, Automatica 
45, 1847 (2009). 

[22] R. L. Hudson and K. R. Parthasarathy, Commun. Math. 
Phys. 93, 301 (1984). 

[23] H. J. Carmichael, Phys. Rev. Lett. 70, 2273 (1993). 

[24] C. W. Gardiner, Phys. Rev. Lett. 70, 2269 (1993). 

[25] A. Barchielli, in Open Quantum Systems III: Recent De- 
velopments, S. Attal, A. Joye and C.-A. Pillet, eds. 
(Springer, Berlin, 2006). 

[26] D. Simon, Optimal State Estimation: Kalman, H x , and 
Nonlinear Approaches (Wiley-Interscience, 2006) . 

[27] J. Gough and M. R. James, Comm. Math. Phys. 287, 
1109 (2009). 

[28] J. Gough and M. R. James, IEEE Trans. Auto. Control 

54, 2530 (2009). 
[29] C. W. Gardiner, A. S. Parkins, and P. Zoller, 

Phys. Rev. A 46, 4363 (1992). 
[30] C. W. Gardiner and P. Zoller, Quantum Noise: A 

Handbook of Markovian and Non-Markovian Quantum 

Stochastic Methods with Applications to Quantum Optics 

(Springer, 2004). 
[31] L. Bouten, R. Van Handel and M. R. James, SIAM J. 

Control Optim. 46, 2199 (2007). 
[32] H. I. Nurdin, M. R. James and A. C. Doherty, SIAM J. 

Control Optim. 48, 2686 (2009). 
[33] F. Marquardt and S. Girvin, Physics 2, 40 (2009). 
[34] T. Botter et al, Phys. Rev. A 85, 013812 (2012). 
[35] M. Tsang and C. M. Caves, Phys. Rev. Lett. 105, 123601 

(2010). 

[36] O. Arcizet, T. Briant, A. Heidmann and M. Pinard, Phys. 
Rev. A 73, 033819 (2006). 



* Electronic address: rhamerly@stanford.edu 



